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1. Abstract 

An Important Issue In the control of practical systems Is the effect of 
model uncertainty on closed loop performance. This Is of particular concern 
when flexible structures are to be controlled, due to the fact that states 
associated with higher frequency vibration modes are truncated In order to 
make the control problem tractable. In — thi s pape r emfr fr ey digital 

simulations of a single-link manipulator system ^ to demonstrate that passive 

damping added to the flexible member reduces adveF'jr^rfectS JSSdc'IA'Teb witn / 

model uncertainty. A controller was designed based on a model Including only 

one flexible mode. This controller was applied to larger order systems to 

evaluate the effects of modal truncation. Simulations using an LQR design 

assuming full state feedback Illustrate the effect of control spillover. 

Simulations of a system using output feedback Illustrate the destabilizing 
effect of observation spillover. The simulations reveal that the system with 
passive damping Is less susceptible to these effects than the untreated case. 



/c 


2. Introduction 


Many In-space robotic operations will require arms capable of very long reach, while like other space 
structures, they must be lightweight. Because such arms are likely to be highly compliant (as Is the space 
shuttle RMS arm), control strategies designed to accommodate structural flexibility must be considered. 
Controlling flexible structures through purely active measures can be cumbersome In terms of hardware and 
computation time requirements. Moreover, active controllers for flexible structures are subject to 
Instability and other problems associated with model uncertainty. The burden of active control can be 
reduced by augmenting active control with passive damping. This enhances system stability and reduces the 
adverse affects of model uncertainty, thereby providing justification for the use of low order dynamic 
models and controllers. 

In this paper we consider a single-link, single-axis arm which rotates in the horizontal plane about a 
pinned hub In response to a control torque t(t) . The system, Illustrated In Figure 1, and the models 

employed In this Investigation are based upon a laboratory version of the arm that has been used In 
experimental Investigations [1-4] at Georgia Tech. The flexible member Is a long slender beam that is 
assumed Infinitely stiff In vertical bending but flexible in horizontal bending. 



307 




Tht pinned hub has rotiry Inertia J. A point payload mass m Is fixed to the bean's tip. Although, 
manipulators sonatinas carry payloads that have significant rotary Inertia, the effect of this Inertia If 
qualitatively similar to that of a point nass for the configuration considered, and hence payload Inertia If 
not Included here. Friction In the pinned Joint Is represented by a rotary viscous dashpot. TMi 

configuration Is viewed as being representative of lightweight, large payload capacity Manipulators, 
Parameters and dimensions for the am that the system considered In this paper are tabulated In Appendix A. 

Damping augmentation Is provided by a constrained viscoelastic layer damping treatment [2,5]. ta 
approach Involves bonding a thin film of viscoelastic material to the flexible member's surface. TMi 
viscoelastic layer In turn has a stiff elastic constraining layer bonded to Its surface. The combine* 
system forms a sandwich-1 ike structure Illustrated In Figure 2. When elastic deflection of the structure 
occurs, shear Induced plastic deformation Is Imposed In the viscoelastic layer. The energy dlsslpatlo* 

associated with the plastic deformation provides the desired mechanical damping. The damping ratio for th* 
untreated beam was approximately constant for all modes at .007. The treatment Increased the damping ratio! 
associated with the modes of Interest (say the first six modes) by about an order of magnitude. The treate* 
beam had a damping ratio of .03 for the first mode and the values for the 2nd through 6th modes ranged fro* 
.052 to .06. Additional damping Is Introduced by joint friction. 


HeetU Cenatr 



Layer 


Figure 2. Treated Beam Clement Under Flexure 


The first step of controller design Is usually the development of a "design model" that Is a slmpllfle 
representation of the actual plant dynamics. The design model serves as the basis for controller design 
In the case of flexible mechanical systems, the design model Is often a truncated representation of th 
actual plant, retaining only a few critical modes. This implies the assumption that a model based upon 
small number of vibration modes provides adequate representation of the much larger order actual plant, fo 
controller design purposes. The modeling error associated with the neglected modes, adversely affect 
closed loop system performance. In this paper, simulation results are presented to illustrate that the 11 
effects associated with modeling error are reduced somewhat through the addition of passive damping to th 
system. 

We consider a multivariable control system, designed according to the steady state linear quadratl 
regulator (LQR) approach. A four state model Including only one flexible mode and the rigid body «od 
represents the design model. We consider the consequences of controlling larger order plants 

representa tlve of the actual system, with a controller derived for the design model. 

The regulator Is formulated to penalize tip position and control effort. Two cases are considered 
The first assumes that full state feedback Is available. The second case uses output feedback of tl 
position (v L ), tip velocity, hub angle (0) and hub angular velocity. The controller designs are kept slnpl 
to facilitate comparisons between the damped and undamped systems. 1 


3* Dynamic Model 

Linear transfer function models for the system of interest were developed based on the assumption ( 
small bending deflections and small hub angles. Transfer function modeling for similar systems has be« 

discussed by several authors [2, 4, 6-8] and we will not repeat the procedure here. Details on development * 
the model employed here may be found In [2], The transfer function poles and zeros used in thl 
Investigation are tabulated In Appendix B. 


1/ Although It possesses some light structural damping and is affected by joint friction, we shall 
designate the untreated arm as "undamped". 
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In the Inurtit of obtaining « tut* ipici realization fro* t h* transftr functions It It convenient to 
work with the partial fraction expansion f ont giv en by Equation 1. Tha w*t art tha nodal frequencies* 

the w*s are tha damped nodal frtquanclts w ✓ 1 - » the C*s art nodal da up In? ratios and n It tha nunbar 

of flexible nodat represented. Damping di i* to Joint friction hat not baan accountod for In that* transfer 
functions but will Introduced later at a fora of feedback* 2 Tha residues k Q and u 0 correspond to the rigid 
body node* 


lc 
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The subscript x on (s) raprasanu the output variable of Interest. For tha present study four transftr 
functions ware required. Equation 2 summarizes these and defines the noutlon used hare. 


hub angular position' 


0<») ‘ 



hub angular rate 

• 

*0(») 

m 

Gj(s) 

beam tip position 


v L (*) 


G vl<*> 

beam tip rate 




. G w. (,) . 


Here s Is tha Laplace operator and T(s), O(s) and ( s ) denote the Laplace transforms of the Input torque* 
hub angle and tip position variables* respectively. 

Figure 3 Is a block diagram representation of the transfer function. 



Figure 3. Block Diagram Equivalent of Equation 1 
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THt tystMi Mtrlcts corresponding to Figure 3 «r« as follows: 


1 

C 

\ o 

** 2 - 2 %<« 

I • [ 0 1 0 1 ••• 0 1 } T (3) 

Si “ C M o x o “I S •" *S» N . 1 

Tho loading principle 2x2 submatrlx of A represents the rigid body node In the tip and hub position transfer 
functions. 

For each of the output variables x there Is a unique output matrix . He will denote thase as 

C e , Cj, C yL and C* L with the noutlon carrying the obvious Meaning. These are assembled to forts a 
Measurement matrix representing four outputs as Indicated In the measurement equation (4) given by: 

* 

C 9 

C i 

, * *C* («) 

l 'yt 

C \ 

In order to account for viscous joint damping we consider the feedback system Illustrated In Figure 4. 





Figure 4. Block Diagram Illustrating Feedback of Joint Damping 


When the effect of joint damping Is Introduced through boundary condition C eedback. a new A matrix Is 
formed: 

A « A - BCjb (5) 


Where b Is the joint damping coefficient. The analysis that follows Is based upon a system model of the 
form (A, B, Q) . 


4. System Representation 

We wish to design a controller for a plant with 2(n+l) states, using a design model Including only one 
vibration mode. Here n represents the number of flexible vibration modes required to provide accurate 
representation of the actual plant, and the additional two states represent rigid body motion. Controllers 
developed for the single mode system are applied to a model Including three vibration modes (8 states), and 
one Including six vibration modes (14 states). These are referred to as the plant models In the text that 
follows because they are Intended to represent actual plants In the simulations presented here. 
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Tta tyttM tQuatlMt for tho flout aodolo art |1*on by: 

i„(t) • ^ (t) ♦ l„« (t) «) 

* (t » * So % (t » 

Here ^ It tut J metrlx defined In Equation 5, la tMt cm for a* a mode approximation of the plant. 
Similarly tl« ^ aatrli Is the ( matrix of Conation 4 for an n mode approx 1 notion. 

THa system equations for tHa design noOals art given by: 

^(t) • ^ «j(t) ♦ IjiKt) (7) 

y(t) • fij ajCt) 

The system output vector (y) considered It of the fom 

y • [ e S v^ v^ ] T . (S) 


S, noftlatar Design 

The standard fort of 1 1noar quadratic cost function Is 


J c ' f (,T °* * " Tr * ) dt 


(») 


where 0 Is a symmetric, posltlvo seml-deflnl te stato weighting natrlx and R Is a symmetric, position 
doflnltt control tffort weighting natrlx. Bocauso no sock to regulate tip position, a porfomanco Index 
that penalizes tho tip position output variable and control effort was chosen. A cost function for tip 
position output weighting Is expressed as follows: 


J 


c 



♦ rt*| dt 


( 10 ) 


Tho output weighted performance Index (10) Is equivalent to the standard stato weighted version (9) with 
weighting matrices given by: 


0 




s« 


R - tr] 


(ID 


The system (A ft , B n , Q n ) represents an actual plant with dynamics that are either incompletely known or too 
cumbersome to permit the use of the full model In controller design. The four stato design 
model (A|, Bj, £j) will serve as an approximation to the actual plant for controller design purposes. In 
this case (j replaces £ n In the state weighting matrix Q (11) 

Two attractive features of the rip position weighted cost function (10) are that It has only one 
parameter (r) to vary, and that a given value of r can be expected to Impose similar performance demands on 
both systems (damped and undamped). Reducing the value of r decreases the penalty on control effort and Is 
therefore equivalent to demanding higher performance at the expense of Increased control energy. 


f. State Feedback 

In a system with decoupled modes, such as the Jordan canonical realization of Equation 3, a state 
feedback law 


u • (12) 

designed to stabilize the reduced order system (7) will stabilize the actual system (6), provided that the 
truncated modes are asymptotical ly stable. The neglected modes can, however, be excited at their natural 
frequencies In response to the applied control Input. This effect Is called control spillover [9,10] In the 
literature related to controlling flexible spacecraft. The system we are considering has a small amount of 
modal coupling, due to the Introduction of viscous joint damping using Equation 5. Since wt normally do not 
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expect viscous damping to destabilize • tystts, It Is reasonable to expect that the state feedback lav (12) 
will stabilize tf» plant models of Intorost. 

In tbls soctlon we design a control ltr based on tba output weighted performance Index (10). We at sum 
that tha first four elements of tha stata vac tor ara somehow available for feedback. Control Is applied to 
these four suta elements according to Equation 12. Tha time varying gain Kj(t) that minimizes the cast 
function (f) is given by 

KjU) - «" ! l 1 T S(t) (13) 

wham S(t) Is tha solution of tha associated matrix Klccatl equation. An often used substitution for tha 
opt leal gain Ki(t) Is Its consunt steady suta value K. • K .(•) . The steady suu value ~f SC t) Is the 
solution of tha algabralc Rlccattl aquation: 

-$(-) - $(-)*, ♦ Aj T S(-) - $(-) •,«' l 8 1 T S(.) *0-0 (14) 


Tha steady stau gain solution is used tha simulations presented hem. 


7. Simulations Kith State Feedback 

Figure S Illustrates tha simulated msponsa of tha design modal (7) to a 4.8 Inch step commend, for 

various values of r. The simulations Indicate that both tha damped and undamped systems are capable of 
almost arbitrarily good performance as r Is decreased. In practice, the limit on the msponsa tine Is 
dictated by the strength of tha beam and tha torque limit of the actuator. A theoretical (assuming fnimite 
beam strength and motor torque capacity) limit on the speed of msponsa Is discussed by Schmitz [6]. This 
limit is related to tha non-minis** phase character of the tip position transfer function. Notice mat the 
wrong way start phenomenon typical to systems with non-minimum zeros Is Indicated In the plots. Schmitz 
Interprets tha theoretical response limit as being roughly equivalent to a pure delay associated with the 
Initial period during which the tip moves In the direction opooslte to the control command. 

When the state feedback law (12) Is applied to the larger order systems (Figures 6 end 7), the 

excitation of the second mode of vibration Is readily apparent when r • 1<T*. In the undamped system 
(Figures 6a and 7a) the oscillation takes mom than two seconds to die out. Thus. In the case of the 

undamped system, we find that designing for higher performance (by reducing r) actually results In slower 
response. The excitation of the second mode also occurs In the damped system (Figure 6b and 7b). however, 
It dies out In about 0.8 seconds. Although the performance of the actual plant Is notes good as tmt of 
the design model (Figure 5b), the simulation Indicates that the response time for r* 10“* Is slightly better 
than the lower levels of demanded performnee (larger values of r) considered. This Is In sharp contrast 
with the results of Figures 6a and 7a for the undamped system. This example clearly Indicates that the 

damped system Is less susceptible to control spillover than the undamped case. 

It should be noted that the peak control torque command, when r • 10“*, Is about 4800 Inch pounds. 
This value Is well above the beam's maximum bending moment capacity (~ 175 In. Ibf. based on yield) end Is 
about 60 times greater than the rated torque capacity (85 In.lbf.) of the experimental system's motor. In 
light of these figures, one might argue that control spillover Is not a realistic concern for the system of 
Interest. The author concedes to the somewhat artlclclal nature of this example, however, further 

consideration of the results adds to their significance. Suppose the Initial step command Is scaled down by 
a factor of 50 to about 0.1 Inches. Because the system model Is linear, we know that the corresponding peak 

torque Is about 100 Inch.lbf. This Is a realistic figure for the system of interest. In Figures 6a and 7a, 

peak tip position oscillation amplitude Is about 1.5 Inches. Scaling this figure down by a factor of 50 
gives 30 thousandths of an Inch - a significant value In the context of robot accuracy. • 


8. Output Feedback 

The simulations presented In the previous section were based on the assumed availability of states. 
Practical control systems must depend upon measured outputs for feedback. Frequently the outputs are 
different entitles than the states. In contrast to systems using state feedback, output feedback systems 
are subject to Instability as a consequence of model reduction [9-11] even when the neglected modes are 
asymptotically stable. This effect Is sometimes called observation spillover. 

In this section we follow the steady state LQR controller design approach employed In the previous 
section, however, we Implement the controller using output feedback. The design model (Aj. Bj, Cj) has four 
states and four outputs and the measurement matrix Cj Is Invertible. This allows us to calculate the state 
Xj of the desing model from the output vector y according to: 


x 


1 


- c i 




y 


( 15 ) 
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a) Undaapad Cast 


b) Oaapad Casa 


Flgura 5. Closad Loop Sup Rasponsa Using Stiti Faadback, Ona Vibration Moda Plant Nodal 









In this cist the output feedback control U-.t Is given by: 


« * Htj (16) 

The output feedback law (16) » when applied to the dtslgn modal Is equivalent to stata feedback (12). When 
applied to controlling tha actual model, the output faadback control law (16), expressed In the form of a 
state faadback law Is given by: 

• ■ -*! S, M S„ *„ W) 

Notice that whan appllad to tha plant modal, tha output faadback law receives Information from tha states 
that wart neglected In design. This Is unwanted Input (spillover), and can be viewed as a form of 
measurement corruption. 


9. Simulations With Output Feedback 

Simulation results obtained using output feedback are presented In Figures 8 through 10. Figures 8a 
and 8b are simulations of the design model using the output feedback law (16). When applied to the design 
model the output feedback law considered Is equivalent to state feedback (Figure 5). This case Is preseated 
here as a basis for comparison. When the low order output feedback law (16) Is applied to the actual system 
models (Figure 9 and 10) we observe that the performance Is limited by the onset of Instability. In the 
undamped system, the first vibration mode Is unstable for r ■ 0.01 and r ■ 0.005. The damped system remains 
stable under the same conditions, however, some first mode oscillatory behavior becomes evident as ve 
attempt to design for higher performance. The damped system Is not Immune to tne Instability experienced by 
the undamped case, however, due to Its more favorable open loop pole placement It 1$ more robust. 

Upon comparison of the six mode and three mode systems, we find that the stable responses of the six 
mode plants do not differ noticeably from those of the three mode plants. On the other hand, the divergence 
rate of the unstable oscillations Is greater In the six mode plant (Figure 10a) than In the three mode pUnt 
(Figure 9a). This Indicates that the presence of the higher, neglected modes (4th, 5th and 6th) do affect 
system stability slightly. 

This example Illustrates that the passive damping treatment considered reduces the flexible system's 
susceptibility to observation spillover Induced Instability. The peak torque commanded at the highest 
performance (when the system Is stable) was about 80 In.lbs.. Indicating that the performance demanded ms 
reasonable for the system under consideration. The example employs perhaps the most simplistic of ill 
possible output feedback schemes. Systems employing state estimators also rely on measured data for 
feedback, and they too are subject to Instability due to modeling error. 


10. Conclusion 

One form of modeling error that Is relevant for control of flexible structures results from Ignoring 
high order vibration modes In the process of deriving a design model. The effects of this type of modeling 
error are mnlfested as control and observation spillover. We have presented simulations of multivariable 
control a particular flexible arm to Illustrate that the addition of passive damping yields a system that Is 
less susceptible to these undesirable effects. 

To some degree these results follow Intuition, In that one naturally expects that Increasing the 
damping terms of a system's eigenvalues will provide a more stable system with Improved performance. The 
results presented are Intended to demonstrate the concept of passive damping on an example that Is 
representative of practical lightweight manipulators. 
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a) Undamped Case 


b) Damped Case 


Figure 7. Closed Loop Step Response Using Output Feedback, Six Vibration Mode .-Plant Model 
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APPENOIX A - System Parameters and Dimensions 


Joint Inertia 
Payload Mass 

Joint Damping Coefficient 
Beam Dimensions 
Material 


J ■ 30.2 1n 2 .lbm. 
m * 0.09 lbm. 
b « 0.10 In.lbf. s z 
48" x 3/4" x 3/16" 
6065-T6 Aluminum 


APPENDIX 8 - Transfer Function Poles and Zeros 


Table B-l Undamped System 


Mode 

System Poles 

Hub Angle T.F. Zeros 

Tip Position T.F. Zeros 

1 

-0.0541 

*37.726 

-0.0149 

*32.1129 

*8.3410 

2 

-0.1284 

*318.3456 

-0.0985 

*314.0768 

*45.0741 

3 

-0.2957 

*342.2446 

-0.2853 

*340.7557 

*111.3047 

4 

-0.5769 

*382.4188 

-0.5719 

*381.6939 

*206.9715 

5 

-0.9633 

*3137.6207 

-0.9603 

*3137.1836 

*332.0743 

6 

-1.4532 

*3207.5965 

-1.4511 

*3207.2946 

*486.6130 

Table B-2 

Damped System 




Mode 

System Poles 

Hub Angle T.F. Poles 

Tip Position T.F. Poles 







1 

-.22197 

*37.1601 

-.0176 

*31-9443 

1 -7.2785, + 7.8830 

2 

-.9124 

*317.4515 

-.6202 

*313.2581 

2 -45.9556, +44.7292 

3 

-2.3169 

*341.6716 

-2.2342 

*340.1610 

3 -97.7161, +113.9306 

4 

-4.8477 

*383.2589 

-4.7987 

*382.5101 

4 -263.5548, +216.0952 

5 

-8.5293 

*3142.4761 

-8.5042 

*3142.0184 

5 -389.2245, +351.1999 

6 

-11.7674 

*3219.5993 

-11.7581 

*3219.2929 

6 -554.2140, +519.0317 


316 







